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1 Introduction 



Nonequilibrium thermodynamical processes are supposed to play a crucial 
role in the physics of the early Universe. Traditionally, for the description 
of these phenomena the theories of Eckart [1] and Landau and Lifshitz [2] 
were used. Due to the work of Miiller [3], Israel [4], Israel and Stewart [5, 
6], Pavon, Jou and Casas- Vazquez [7], Hiscock and Lindblom [8] it became 
clear however, that the Eckart type theories suffer from serious drawbacks 
concerning causality and stability. These difficulties could be traced back to 
their restriction to first order deviations from equilibrium. If one includes 
higher order deviations as well, the corresponding problems disappear. By 
now it is generally agreed that any analysis of dissipative phenomena in rela- 
tivity should be based on the theories of Miiller, Israel and Stewart, including 
at least second order deviations from equilibrium, although in specific cases 
the latter might reproduce results of the Eckart theory [9]. Cosmological 
implications of second order theories were first considered by Bclinskii et al. 
[10]. In the realm of cosmology especially bulk viscous phenomena have at- 
tracted considerable interest (see, e.g., [11]) since bulk viscosity is the only 
possible dissipative mechanism in homogeneous and isotropic spacctimcs. 
While the coefficient of bulk viscosity vanishes both for pure relativistic and 
pure nonrelativistic equations of state, it may be important, e.g., for mix- 
tures of radiation and matter [12]. On the other hand, it is well known 
[13-15] and widely used [17-25] that particle production processes in the ex- 
panding Universe may be phenomenologically described in terms of effective 
viscous pressures. 

A major point of interest in the study of bulk viscous universes has 
been the question whether there are conditions under which a sufficiently 
large bulk viscous pressures could lead to an inflationary behaviour. While 
some authors concluded that a bulk viscosity driven inflation is impossible 
[26] , others [27-29] found inflationary solutions of the cosmological evolution 
equations. Partially, these difl'erences occur since different equations of state 
were used. But the results also depend on whether the investigations were 
performed within the full, causal second order theory or in a truncated 
version of the latter. 

One should be aware that discussing the issue of bulk viscous driven 
inflationary solutions at all, implies in any case an extrapolation of nonequi- 
librium thermodynamical theories beyond the range for which their appli- 
cability was strictly justified [29]. Bulk viscous inflation, if it exists, is a 
far-from-equilibrium phenomenon, while even the full causal second order 
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MIS theory is a theory for smah deviations from equihbrium. Therefore all 
theoretical conclusions are necessarily tentative. It is the hope that they 
nevertheless will provide an indication of the correct behaviour far from 
equilibrium. 

As to the relation between the full theory and the truncated version, to be 
discussed in some detail below, the following comment should be made from 
the outset. To decide whether a theory is truncated or not on apparently 
obvious formal grounds, i.e., from the appearance of the causal evolution 
equations may be misleading. The structure of the evolution equation de- 
pends on the choice of the basic thermodynamical variables. In most cases 
the latter are equilibrium variables and the Gibbs equation has their fa- 
miliar form. In the framework of 'Extended Irreversible Thermodynamics' 
(EIT) however, a generalized Gibbs equation is used which includes dissipa- 
tive quantities as independent variables. Temperature and pressure in EIT 
are nonequilibrium quantities, different from their equilibrium counterparts. 
Written in terms of these nonequilibrium quantities the causal evolution 
equation formally may look identical to the truncated theory written in 
terms of the more familiar equilibrium variables. These points have been 
clarified in a recent paper by Gariel and Le Denmat [30] who pointed out 
that the apparently (because of the use of nonequilibrium variables) trun- 
cated theory of Pavon et al. [7] in fact is equivalent to the full theory of 
Israel [4] . The reader is also referred to corresponding comments in [28] and 
[29]. 

In the present paper the symbols T and p always denote the equilibrium 
temperature and the equilibrium pressure, respectively. 

In section 2 we reconsider the bulk viscous cosmological dynamics within 
the causal second order theories. The evolution law for the temperature of 
the cosmic fluid is shown to be different in general from the ad hoc relation- 
ship used in previous treatments by Romano and Pavon [31], Zakari and Jou 
[28] and Maartens [29] (subsection 2.1). In subsection 2.2 we discuss the con- 
ditions under which the truncated version yields results close or identical to 
those of the full theory. The general viscous fluid dynamics of the full MIS 
theory is presented in subsection 2.3. Subsection 2.4 investigates the con- 
ditions for viscous exponential inflation. Almost all corresponding solutions 
imply an exponential increase of the fluid temperature. As a consequence 
of this behaviour there is in general no substantial growth of the comov- 
ing entropy as was found previously [29]. The latter result corresponds to 
a very specific limiting case. In section 3 the viscous pressure of the full 
causal theory is assumed to describe partially or fully the effect of particle 
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creation taken into account by a nonvanishing source term in the particle 
number balance. In this setting the backreaction of the viscous pressure 
on the temperature is different from the conventional viscous fluid case of 
section 2. In the limit that the viscous pressure is entirely due to particle 
production there exist stable inflationary solutions for which both the par- 
ticle number density and the temperature are constant and, moreover, the 
comoving entropy grows exponentially. Section 4 summarizes the results of 
the paper. Units have been chosen so that c = fcs = 1. 



2 Bulk viscous fluid dynamics 
2.1 General relations 

The energy momentum tensor of a relativistic fluid with bulk viscosity as 
the only dissipative phenomenon is 

T'^ = pu'u^ + (p + tt) . (1) 

p is the energy density, is the 4- velocity, p is the equilibrium pressure, 
is the projection tensor = g^^ + u^u^^ and tt is the bulk viscous pressure. 
The particle flow vector iV" is given by 

AT" = nn" , (2) 

where n is the particle number density. Limiting ourselves to second order 
deviations from equilibrium, the entropy flow vector S"" takes the form [5, 
26] 

2 

2CT ^ ' 

s is the entropy per particle, r is the relaxation time, T is the temperature 
and C is the coefficient of bulk viscosity. The conservation laws 



and 

imply 

and 



n + en = , 

p = -e (p + p + tt) , 



(4) 

(5) 
(6) 

(7) 



4 



respectively, where = is the fluid expansion and n = n^au"" etc. Com- 
bining (|6|) and (0) with the Gibbs relation 



Tds = d— + pd— 
n n 



we get 

From (^) and (|6|) we find 
S^. = - 



nTs 



-en 



e + -7T + -vrT — 



C 



CT 



(8) 
(9) 

(10) 



for the entropy production density 5.^. The simplest way to guarantee 
5." > implies the evolution equation 



vr + Tvr 



e + 



f c 



r C 



T 
T 



for vr, leading to 



(11) 



(12) 



For r — > eq.([ll|) reduces to the corresponding relation of the Eckart theory. 
The frequently used truncated version 



vr + rvr 



-ce. 



(13) 



also known as Maxwell-Cattaneo equation, follows if the bracket term on 
the r.h.s of (|ll|) can be neglected compared with the viscosity term — C©- 
Below we shall give explicit criteria for this approximation. If the bracket 
term vanishes identically, i.e., if the condition 



e + 



r C T 







(14) 



is fulfilled, the full and the truncated theories become identical. As we shall 
see this is possible only in exceptional cases. 

Let us assume equations of state in the general form 



p = p{n, T) 



(15) 
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and 

p = p{n,T) , (16) 

according to which the particle number density n and the temperature T 
are our basic thermodynamical variables. Differentiating the latter relation, 
using the balances (^) and (0) as well as the general thermodynamic relation 

dp _ p + p T dp ^^^^ 
dn n n dT ' 

one finds the following evolution law for the temperature (cf. [22]): 

■ dp/dn p 



or 



T 



dp/dT ^ vr 



dp/dT Tdp/dT 



(19) 



For TT = and with = 3R/R, where R is the scale factor of the Robertson- 
Walker metric, (|l^ ) reproduces the well known Tj. ~ R~^ behaviour in a 
radiation dominated Friedmann-Lemaitre-Robertson- Walker (FLRW) Uni- 
verse, while for p = nm + |reT, p = nT and T <^ m we recover T^ ~ R~^ 
in the matter dominated case. For a viscous fluid the behaviour of the tem- 
perature depends on vr. Since vr is expected to be negative, the second term 



in the bracket on the r.h.s. of (19) will counteract the first one. Close to 
equilibrium, i.e., for | vr |< p the existence of a bulk viscous pressure implies 
that in an expanding universe the temperature decreases less rapidly than 
in the perfect fluid case. 



2.2 The truncated version 

While the truncated version was used in most of the earlier applications, 
more recently an increasing number of authors [26, 28, 31, 29] has studied 
the full theory and compared the results of the latter wih those of the trun- 
cated version. In some cases these results differ dramatically, which may 
be interpreted as a breakdown of the Maxwell-Cattaneo type equations as 
a reasonable approximation to the full theory under the corresponding con- 
ditions. What seems to be missing, however, are general criteria according 
to which one may decide whether the truncated version is sensible and be- 
yond which limits it fails to give an answer close to that of the full theory. 
Intuitively one expects the coincidence to be the better the closer one is to 
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the equilibrium case. Below we shall give an example showing that there 
are identical results in exceptional cases even far from equilibrium. In order 
to clarify the approximative character of the truncated theory we assume, 
as usual, the relation C, = pr that guarantees a finite propagation velocity 
of viscous pulses [10, 27-29, 31]. In the following subsection we are going to 
generalize this relation. Using (0) and (19) in equation (0) with ( = pr we 
find 



vr + rvr = — pGr 
with the abbreviations 



l + ^(l + 6 + 7) + ^(a + l) 



(20) 



1 dp ^_dp/dT I dp p + p 



nTc^ ndT dp/dT nc^ dT p 

Obviously, the truncated version is expected to be applicable for 



vr I , . , vr^ 



1+6 + 7) < 1 , 77^(1 + 1 • (22) 



2p ^ " ' 2p 

Since for 'ordinary' matter h and 7 lie in the ranges 1/3 < b < 2/3 and 
1 < 7 < 4/3, respectively, the first condition is roughly equivalent to it <^ p. 
For radiation with 7 = 4/3, b = 1/3 and a = 1 the second condition is 
implied by the first one. For matter with 7 = 1,6 = 2/3 and a = |^ ^ 1 
however, the second condition has to be checked separately. The inequalities 
(p2| ) may be regarded as criteria under which the truncated theory is a 
reasonable approximation to the full theory. Given equations of state (15) 



and ([l6|), any solution vr of the truncated theory may be tested according to 
(22) whether or not and to which accuracy it approximates the full theory. 
But there is a different possibility, namely the case 



2 

vr , . , vr 



-(l + 6 + 7) + ^(a + l) = 0, (23) 

in which all the terms that distinguish the full from the truncated theory 
cancel among themselves. Relation ( p3| ) is identical to ( p^ for ( = pr, i.e., 
it is the condition under which the full theory is identical to the truncated 
one. Solving ( p3[ ) for vr//? yields 

!i = _l±2±i:. (24) 

p 1 + a 
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According to the above mentioned parameter ranges for 7, b and a, | tt |<C p 
is only possible for a 3> 1, i.e., for massive particles. Using (p^) in (|^ yields 

^ = -e^^^^^i^ . (25) 

p l+a ^ ' 

For radiation (7 = 4/3, b = 1/3, a = 1) we find p = and, consequently, 
G = 00 = const in a flat FLRW universe. The truncated and the full theory 
coincide in a specific bulk viscosity driven inflationary universe. Since with 
vr is completely determined by p, provided, the equations of state are 



given, the remaining equation (13) is no longer a dynamical equation on 
its own, but may be used to calculate r = r(Go). Since the solution is 
stationary, i.e., vr = 0, we find 

rS = l^o' ■ (26) 

In a radiation dominated universe the full theory and the truncated version 
admit a common, bulk viscosity driven inflationary solution with a relax- 
ation time of the order of the expansion time. This seems to be a new result. 
For the temperature dependence we find from (^) and ( p^ 

T l+a ' ^ ^ 

yielding 

^ = e (28) 

for radiation, or T ~ i?'^. The temperature increases in an expanding uni- 
verse. This implication of the condition (0) was first noticed by Maartens 
[29]. While this might appear strange at the first glance, it is an unavoidable 
feature of any bulk viscosity driven inflation as is obvious from (|l9|). This 
point will be discussed in more detail in the following subsection. 

2.3 The dynamics of the full second order theory 

In this subsection we shall investigate the full causal theory assuming the 
existence of general equations of state ([T5|), (p^) and 

-^f = fip) ■ (29) 
r 
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(30) 



Following Belinskii et al.[10] usually the relation C/t = p was used to guar- 
antee that the propagation velocity of viscous pulses, which is expected to 
be of the order [28] 

does not exceed the velocity of light. We shall not immediately specify to 
/ = /9 in order to admit a certain range for this propagation velocity. As 
we shall see below, this additional freedom, allowing, e.g., f = ap with 
< a < 1, may be useful in characterizing a possible inflationary phase. In 
this case one has / = f'p, where /' = df /dp. Consequently, using 



f_ 

f 



TT 



@Pj ( 7 + - 



Together with (^), equation (jTT|) may now be written as 



TT + TTT 



-pQt 



f 



IT 



p 2p 



Pf 



- + — 1 + + 7^ + TTT a + ^ 



/ 



2p^ 



pf 



f 



(31) 



(32) 



It is obvious how the applicability conditions (|22| ) of the truncated version 
have to be modified in this more general case. Restricting ourselves to a flat 
FLRW universe with 

02 



where k is Einstein's gravitational constant, and 

3k 



G 



(p + P + tt) , 



(33) 



(34) 



the latter equation may be used to eliminate vr and vr, respectively. The 
bulk pressure tt may be written as 



-3jH^ - 2H 



where 3H = Q and H = R/R or. 



TT 
P 



-7 



2 H 
3lP ■ 



From (H), using (§) and (|T9|) one realizes 

P = vIp + Qtt [vl - bj 



(35) 



(36) 



(37) 
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with the sound velocity Vs, given by 
'dp 



n 



dp T {dp/dry 



Differentiating 
replaced by 



KTT 



dp/ isentropic P + P On p+p dp/OT 
and using and (|37|), ktt on the l.h.s. of (| 

= -2H - 6HH (1 + 6) + 9H^-f (v^, - 



(38) 
may be 

(39) 



Applying 
for H: 

tH - 



and (By) in (p3) we arrive at the follwing evolution equation 



El 

H 



3HHt 



, 1 P fi 
7a H f 



2 



- + ^ (7« - 1 - f*) + 7 U 
p 2 V 



;(l + 6) 

2 ' 



. (40) 



It should be pointed out again that in arriving at this evolution equation 



only equations of state (|l5D , (|lq) and (|29|) were used. We think this set of 
equations of state to be more general than those used in previous papers 
[31, 28, 29]. Especially, no specific temperature law like T = Pp^' in [31] had 
to be postulated. There is no freedom to impose a separate temperature 
law. The behaviour of the temperature is generally governed by (19). As it 
is obvious from (|l^), a relation T ~ p^' is possible for dp/dn = 0. Even in 
this case, however, r is not arbitrary but determined by r = p (Tdp/dT)^^ . 
In the specific case of a flat FRLW universe with 



H 
H 



2p 



equation (^) reduces to 



- = 3(7a-6)- + a-. 

For constant values of a, 6 and 7 

T ~ p"" R^('r'i-b) 

results, which in the radiation dominated case with a 
specifies to 



(41) 



(42) 



1, b 



(43) 

1/3, 7 = 4/3, 



Obviously, with pr 



Tr ~ PrR-^ . (44) 

R~'^ for TT = 0, one obtains the correct limiting case 
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2.4 Viscous exponential inflation 

It has been a matter of some debate whether a bulk viscous pressure is able 
to drive inflation [26-29, 31]. While Hiscock and Salmonson [26] have shown 
that, using the equations of state for a Boltzmann gas, the full causal theory 
does not admit an inflationary phase, different equations of state like those 
used by Romano and Pavon [31], Zakari and Jou [28] and Maartens [29] are 
compatible with a de Sitter phase. It should be stressed again that dealing 
with the question of a viscosity driven inflation one has to assume that the 
MIS theory is applicable far from equilibrium [29]. Under these premises 
we are going to consider now the possibility of a solution of the evolution 
equation for H with H = Hq = const. In such a case eq.(^0|) yields 



— + Uia-l-b)+vl 
IP 2 



(45) 



In the radiation dominated case this reduces to 

r'o = IjH^' ■ (46) 

Keeping in mind that until now only f = f {p) was used, one recognizes that 
the ratio between Tq and Hq^ crucially depends on f / p. If, as is usually 
assumed, f = p, we have Tq < H^^. If, however, / is allowed to be smaller 
than p, we may have Tq > Hq^. It follows from (|30| ) that a relation f < p 
implies a lower propagation velocity compared with the case f = p- Provided 
f = p is equivalent to a propagation velocity that coincides with the velocity 
of light, f = p/3, e.g., leads to a propagation with Vs = l/\/3, the velocity 
of sound (see Israel and Stewart [6]). In the latter case the relation ( |46| ) 
specifies to 

r'o = ^H^' , (47) 

implying Tq > Hq^. In this case the nonequilibrium is 'frozen in'. The 
viscous inflation is 'nonthermalizing' [29]. While for Tq < H^^ there is a 
quick (compared with the expansion rate) relaxation to (local) equilibrium, 
this is no longer the case if Tq becomes comparable to the expansion rate 
Hq . Nonequilibrium situations like this may occur, e.g., in GUTs close to 
the Planck time where the underlying microscopic process is the decay of 
heavy vector bosons [32]. According to Maartens [29] the condition tq > 
Hq^ may be regarded as a consistency criterion for causal viscous inflation. 
Clearly, for tq < , the Universe will relax to an equilibrium state in less 
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than one expansion time. A successful inflation, however, has to last for 
many characteristic expansion times. 

In the matter dominated case with 7 = 1, 6 = 2/3, ?7^^1,a = |^^>l, 
we find 

~ ^H,' , (48) 

i.e., we are always in the range Tq" ^ . The time scale for the relaxation 
to equilibrium is much smaller than the expansion time scale. 

Calculating the dependence tq (-f^o) along the the same lines that lead 
to (]45| ) within the truncated version provides us with 



37-^ 



/ , f 2 

O V 



Ho . (49) 



Specializing to the radiation dominated case, (l49|) reduces to (l4q), i.e., the 
truncated version yields exactly the same dependence as the full theory. We 
have rediscovered the specific case, discussed in subsection 2.2, in which the 
full theory coincides with the truncated version. For the matter dominated 
case, however, the result of the truncated theory is different from (|4^). This 
is not unexpected as well after the considerations of subsection 2.2. 

We shall now look at the behaviour of the basic thermodynamic quan- 
tities under the condition H = Hq. Because of ( ]40|) this implies p = po = 
const and ( [3^ ) reduces to 

- = -7 . (50) 
P 

According to (0) the particle number density decreases exponentially. The 



temperature evolution law (|T9| ) together with (50) becomes 

'dp/dT p 



T 

7f = -3^0 



dp/dT nTcy 



7 



(51) 



equivalent to 



t = ^|^. (52) 
c„ on 



Generally, the temperature is not constant in the inflationary phase. Pro- 
vided c„ is finite, a constant temperature T = Tq is only possible for 
dp/dn = 0, corresponding to T ~ p'' with r = p {Tdp/dT)~' (see the dis- 
cussion below (^)). For any dp/dn > the temperature T increases during 
inflation. Since dp/dn > and dp/dn = 0, respectively, lead to dramatically 
different results for the entropy production during a bulk viscosity driven 
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inflationary phase (see below), this point deserves a detailed discussion. It 
is frequently used that in the equations of state for radiation in equilibrium, 
p = riT and p = 3nT, which are specific cases of ( p^ ) and (p^), the number 
density n may be eliminated according to n ~ T^, resulting in p = p{T) 
and p = p(T) with p = p/3 and p ^ T"^. Then T is the only independent 
variable and dp/dT dp/dT. In order to check whether a corresponding 
procedure is possible for vr 7^ 0, we shall assume an arbitrary dependence 
n = n (T) instead of the equilibrium relation n ^ T"^. With ( [T9D we find 

n = -3HT^(b + a-) . (53) 



dT \ p 

Comparison with (^ leads to the general relation 

dT /, 7r\ dn , ^, 

= 6 + a- — . (54) 
TV pj n 

For vr = we recover T ~ n^^^ for radiation (b = 1/3). For matter (b = 2/3) 
the correct result T ~ i? ^ is obtained as well. 

With the condition (|50|) for bulk viscosity driven exponential inflation 
and for a = const equation (^3) yields 



T ~ n'"-^" . (55) 
In the radiation dominated case 

n ~ T"^ (56) 

results. While in equilibrium n is an increasing function of T, n decreases 
with T in a de Sitter phase, characterized by ([50|), which is a far-from- 
equilibrium state. According to ( |56|) the exponential decrease of the particle 
number density following from (|6|) neccessarily is accompanied by a corre- 
sponding increase in the temperature. This is obviously incompatible with 
T ~ p^/^ from dp/dn = 0. Equations of state with dp/dn = imply that 
the exponential dilution of the particles of the out-of-equilibrium fluid ac- 
cording to (^) does not have any impact at all on the energy density of this 
fluid. We conclude that the previously [31, 28, 29] used ad hoc assumption 



r ~ p'' is not consistent with (p^), (16) and n = n (T) if applied to a bulk 
pressure driven inflationary phase. 

While an increasing temperature during the de Sitter stage appears un- 
familiar, this kind of behaviour is not quite unexpected for an equation of 
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state (p!6[). Since n decreases exponentially, T must increase accordingly in 
order to guarantee p = 0, as long as dp/dn > and dp/dT > 0. Exactly this 
kind of behaviour was found by Hiscock and Salmonson [26] for the trun- 
cated Israel-Stewart theory. As was demonstrated above, the full theory 
exhibits a corresponding feature as well. (Note, that according to [26] there 
does not exist a de Sitter phase in the full theory for the case of a Boltzmann 
gas). If one does not use, however, a relationship T ~ p'' and the behaviour 
of T is governed by ( |52D with dp/dn > 0, this has important consequences 
for the entropy production during bulk viscous driven inflation, discussed 
by Maartens [29]. The entropy in a comoving volume is S = nsB?. With 
(^) and (P) the change of S is 

t = -^-^R'. (57) 
In the de Sitter phase with (|^) the latter expression reduces to 

tH=H, = ^Hoj!^ . (58) 

This change of the comoving entropy depends crucially on the behaviour of 
T. For T = const which follows from T ~ p'^ , equivalent to dp/dn = (see 
(|52|)), we recover the exponential increase of S found by Maartens [29]. For 
equations of state (|l5|) and (|l6|) however, with dp/dn > 0, the temperature 
dependence in the case of radiation (7 = 4/3) is Tr ~ PoR^ as follows 
from (^). S turns out to be constant rather than exponentially increasing. 
Integrating ( [58| ) one finds only a linear growth in S. Consequently, there is 
no way to generate a considerable amount of entropy during a bulk viscosity 
driven de Sitter phase. 

To investigate the stability of the solutions H = Hq, we probe the latter 
with small perturbations, i.e., we assume H = Hq {1 + h(t)) with | /i [^ 1 in 
eq.(|^). Equation (pO|) is valid for spatially flat homogeneous and isotropic 
spacetimes, therefore the stability analysis is restricted to this case as well. 
We assume that the dimensionless quantities f/p, 7, a, b and Vs remain 
unchanged for small deviations from H = Hq. Since we found the relaxation 
time To in the inflationary phase to be proportional to Hq^ it is natural to 
assume r ~ H~^ generally and to fix the proportionality factor by (|45|). 
Retaining only terms linear in h, the resulting equation for h is 

h + SHQKh = , (59) 
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/ , ,2 '1 f P 



with 

IP 2 V/ / 

There exists a solution h = const that may be used to redefine Hq. The 
other solution is stable for K > Q. In the radiation dominated case 7 = 4/3, 
v1 = 1/3, a = 1, 6 = 1/3 and with f = ap the solution is stable for any 
a > 0. In the matter dominated case 7 = 1, <C 1, a ^ 1 and b = 2/3 we 
find K ~ —a/2 = — m/3T, i.e., the solution is unstable. 

We conclude that there exist stable viscosity driven inflationary solutions 
as long as the equation of state is close to that for relativistic particles. There 
do not exist stable solutions for equations of state close to that for dust. Even 
for the stable solutions, however, the temperature increases exponentially 
during the de Sitter stage. Since the particle number density decreases 
exponentially, this unfamiliar behaviour is unavoidable to guarantee p = po 
as long as dp/dn and dp/dT are assumed to be positive. A temperature 
dependence like this is not as strange as it might appear at the first glance. A 
nonvanishing bulk viscosity has always the tendency to heat up the Universe. 
Close to local equilibrium this means that the decrease of T due to the 
expansion of the Universe is less than without bulk viscosity. Applying the 
cosmological dynamics with bulk viscosity to situations far from equilibrium, 
in our case to a hypothetical inflationary phase, the vr-term just compensates 
the equilibrium terms in (0) and (|3j). It overcompensates them, however, 
in the case of the temperature (19). 

In the following section we show how particle production processes may 
modify this behaviour of the thermodynamic quantities during a de Sitter 
phase. 



3 Bulk viscous pressure and particle production 

3.1 Basic dynamics 

Throughout this paper the Universe is studied within a single fluid model. 
There is entropy production due to a nonvanishing bulk pressure. While 
both for pure radiation and pure dust the bulk viscosities tend to zero, con- 
siderable values of the latter are expected in mixtures of relativistic and 
nonrelativistic matter. Consequently, the single fluid universe with bulk 
viscosity may be regarded as a simplified description of a system of two 
(or more) interacting components. Usually it is assumed that the interac- 
tion responsible for the existence of a nonvanishing viscous pressure does 
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not change the overall number of fluid particles, i.e., it is assumed that (^) 
holds. However, interactions with conserved particle numbers are only a spe- 
cial case, particularly at high energies. Moreover, at times of the order of the 
Planck time or during the scalar field decay in inflationary scenarios, par- 
ticle or string production processes are supposed to affect the cosmological 
dynamics, leading to features like 'deflationary universes' and 'string-driven 
inflation' [16 - 20, 38]. 

In this section we generalize the previous formalism to the case that the 
existence of a nonvanishing bulk pressure is accompanied by an increase in 
the number of fluid particles. We start our investigations with (||), (^), (^) 
and (H) as well, but @ is now replaced by 



yielding 



Nt. = nV 



h + Qn = nT . 



(61) 



(62) 



r is the particle production rate which has to be regarded as an input 
quantity in our phenomenological description. Instead of (^) now 



nTs = —Qtt — {p + p)T 
results. The entropy production density is given by 



T 1 

C 2 \CT 



where is the chemical potential 



ji = Ts 



P + P 



n 



(63) 



(64) 



(65) 



A change in the number of particles, i.e., F 7^ 0, is believed to be phe- 
nomenologically equivalent to an effective viscous pressure [13, 15, 17 - 21]. 
A discussion of this equivalence on the level of relativistic kinetic theory 
has been given recently [33]. In order to relate the present investigations to 
previous work we shall in a first step focus on the Eckart theory, i.e., to the 
case r = 0. Equation (p3) reduces to 



E\a 



-npT — tteQ 



(66) 



in this case, where the subscript E stands for 'Eckart'. Formal rewriting 
yields 



+ rifiTT^^T 



(67) 
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As usual, one has to guarantee S^.^^ > 0. Since we expect the particle 
production to be effectively equivalent to an viscous pressure we have to 
demand additionally that the entire entropy production is given in terms of 
TTE, i.e.. 



oa 



where tte is determined by 



< 



+ nfiir^^r 



(68) 



(69) 



or, 



(70) 

where (" now is a generalized bulk viscosity coefficient. For T ^ the latter 
inhomogeneous quadratic equation replaces the familiar linear relation tte = 



— C©. Via (70) the particle production rate T influences the viscous pressure 
TT. For r = as well as for = we recover tte = — C©. It is convenient to 
split in the following way. Let A be a not necessarily constant parameter 
lying in the range < A < 1 such that the fraction Xtte of tte describes a 
'creation' pressure, due to F 7^ 0, while the fraction (1 — A) tte is connected 
with a conventional bulk viscosity. With this splitting (|6^) may be written 



as 



riTsE 



[1-X) Otte- XQtte- {p + p)T . 



(71) 



Frequently [22 - 25] the assumption was made that the creation process does 
not affect the entropy per particle. This is equivalent to the requirement 
that the terms in (71) due to the creation process cancel among themselves: 



{p + p)T = -XtteQ . 



(72) 



Physically this means that the particles are created with a fixed given en- 
tropy. Then S£; is given by 



nTsE = -(1 - A) Otte 



(73) 



Only that fraction of vr^; that is not connected with a change in the particle 
number contributes to s^;. In the limit A = 1, in which the entire viscous 
pressure is due to particle production, one has se = and S'j^.^ = ns£;F, i.e., 
there is entropy production because of the enlargement of the phase space. 
Using (^) in (|6|) yields 



TTE 



-CG 



1 - A- 



P + P. 



(74) 
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Combining the latter relation again with ( [72D we find that the part of the 
generalized coefficient of bulk viscosity C that arises due to a nonvanishing 
particle production rate is 

re-2 

AC = ^ , . (75) 

Of course, ( [7^ ) with (|^) is also obtained if the condition ( |72[ ) is immediately 
introduced in (]66|). 

Obviously, (|6^) with (|6^) is not the only way to guarantee 5." > 0. For 
a different approach that treats F as a thermodynamic flux independent of 
TT, the reader is referred to a paper by Gariel and Le Denmat [34]. 

Coming back now to the full second order theory again, we try to find 
a causal evolution equation for vr that yields (|T^ ) in the case F 7^ as well. 
Equation (E3) may be written as 



W = -TT 



e + -vi- + -ttT —u" + ^ — 



(76) 



C 2 VCT ),a ^ 
generalizing the relation (^) of the Eckart theory. The condition > 



with (12) is now fulfilled for 



vr^ + rvr^ + ]-ti'^T (^u''\ + CvrO = -C//nF , (77) 

instead of (|70|). The viscous pressure vr is determined by a nonlinear inho- 
mogeneous differential equation. For /i = eq.(|7^) formally reduces to ( |lT| ) 
For the sake of generality we shall however retain the //-term in the following 
considerations. A chemical potential may act, e.g., as an effective symmetry 
breaking parameter in relativistic field theories [35-37]. 

In the Eckart theory the general nonlinear relation (^0|) was reduced to 
the linear relation ( |7^ by the requirement (£2). Unfortunately, a corre- 
sponding simplification of ( [77|) via a relation (|72) is not possible in general. 
The essential physical difference between the noncausal and the causal the- 
ory is the appearance of a finite relaxation time within the latter. If a 
nonvanishing F is responsible for the occurence of an effective viscous pres- 
sure, a causal theory has to include a finite relaxation time which is just the 
time interval during which the corresponding part of vr decays to zero after 
F has been switched off. But the relation (^) is an Eckart type relation due 
to which Avr = is implied by F = immediately. Therefore an approxi- 
mation (^) cannot be used in a causal theory. If however, the relation ([7^) 
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which fohows from the requirement that the creation process does not affect 
the entropy per particle has to be abandoned, this means that within the 
causal theory the creation process contributes to s in general. While in the 
Eckart theory the particle production rate was rather simply related to the 
viscous pressure by ( [7^ ) , F enters the equation for vr of the causal theory as 
an independent parameter in a less obvious way. 

Our present approach assumes that the deviations from equilibrium may 
be characterized in terms of one single quantity vr also for the case of a non- 
vanishing particle production rate T. This implies that there is only one 
(generalized) coefficient of bulk viscosity C and one relaxation time r. The 
advantage of this assumption is that there exists only one causal evolution 
equation, namely (^), taking into account the influence of T on vr. (For 
an alternative proposal see again [31]). Since the latter equation is nonlin- 
ear there is no obvious separation into a conventional bulk pressure and a 
'creation' pressure. We shall assume furtheron, that ( and r continue to 
be related by (|29|). The dependence of the relaxation time r on F will be 
discussed below for specific cases. 

For = eq.(|77|) appears to be identical to the case F = 0, but there 
are differences in the behaviour of the particle number density and the tem- 
perature. The former is determined by (|6^). Using (|T6|), ([T7|) , (0) and (|6^) 
instead of (^), the temperature law ([l9|) is replaced by (cf. [22]) 



T 
T 



-e 



dp/dT 

dp/dT ' Tdp/dT 



+ 



+ F 



dp/dT 



p + p 



dp/dT Tdp/dT 



(78) 



in the case of a nonvanishing particle production rate. The particle produc- 
tion affects the temperature not only through the effective viscous pressure 



TT but there is an additional direct coupling as well. Alternatively, (78) may 
be written as 



T 
T 



-(e 



dp/dT iTe + {p + p)T 



(©-r; 



dp/dT 
dp/dT 



Tdp/dT 



+ 



ns 



(79) 



dp/dT dp/dT 

With this evolution law for the temperature, different from (^) the evolution 
equation ( [77| ) for vr becomes different from (32) even in the case /U = 0, 
although (^) and (|7^) seem to coincide. The relation ( p7D for p is now 
replaced by 



p = p 



{T-Q)^v'l-h 



7F + e- 
p 



(80) 
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Instead of f|39l) we have 



-2H - 6HH (1 + 6) + 9H^-/ (w^ - ft) - SH^F-/ (vf _ 5) 



The evolution equation for vr may be written 
TT + rvr = —pr 



p \ p + pir/p 



G (1 + 6 + 7^ 
Pf 



TT 

+ ( + ^ 



2p 



f 



T{b- -fa) 



(81) 



(82) 



generahzing ([34 ) . A procedure analogous to that leading to ( [40D in section 
2 provides us with an evolution equation for H in a flat, homogeneous and 
isotropic universe with particle production: 



3HHt [7a + 1^/'-^ (1 + 6) 

+H + iriJr (7a - b) 
+ l{ja-l-b)+j[vl-b) 



f 
— b 



+ 



37772 



f np 



2 H 1 





(83) 



For an easier comparison with eq.(^) we kept separate the common factor 
ir/p (see (^)) on the l.h.s. of this equation. Having only one evolution 
equation for H isa consequence of our previous assumption that it is possible 
to characterize the deviations from equilibrium by a single quantity vr only. 
Eq.(|83|) simplifies considerably for p = 0: 



3HHt 



2 



+H + irifr (7a - b) 
L + l{^a-l-b)+^(vl-b) 
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+-tH Tj {v,-b + 



3 

2' 



. 



(84) 



For r = both equations reduce to (HO). Eliminating vr in the evolution 
equation (^) for T with the help of (^6|), we find 



^ = 3 (7a - 6) ^ - (7a - 6) r + 
1 It p 



instead of (^21). For constant values of o, b and 7 
For 7 = 4/3, 6 = 1/3, a = 1, this specifies to 



N n 



(85) 



(86) 



(87) 



3.2 Inflationary solutions 



Looking for solutions H = Hq = const of ( 82 ) , we find the following expres- 
sion for To that generalizes (^) 



3Fn 



/ ^ 7Q 
7P 

n/i / 



1 



+ 



2 

7a ■ 



-P + PIP 

If the particle production rate is comparable with the expansion rate, it 
may essentially influence the relaxation time tq. Especially interesting is the 
possibility that the resulting effect of the F-terms in (^) is to enlarge tq. 
In this case particle production may either enable or improve the fulfilment 
of the 'freezing in' condition tq > -f^o^^ for the nonequilibrium. In other 
words, the consistency criterion for inflation might be easier to fulfill in the 
case with particle production than without. Of course, these considerations 
make only sense as long as tq remains finite. As we shall see below there are 
parameter combinations for which tq diverges. 

If during some time interval the particle production rate F is proportional 
to the expansion rate and approximately constant as well, i.e., F = 3XHq, 
we get 



3Ha 



^11 

IP 



Xnp \ 
P + pJ 



A) 



7a 



b 2 



(89) 
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For particles with fi = and an equation of state close to that for radiation it 
is obvious that for any A > the relaxation time tq is larger than for A = 0. 
In the limiting case A = 1 the production rate of the particles coincides with 
the expansion rate, tq then is given by 



1 {A=l) 

-IP \ P + PJ 2. 
Specifying again to 7 = 4/3 and restricting ourselves to /i = 0, 



(90) 



(A=i) 

p 2 

results. For f = p, the usual choice in the literature [10, 28, 29, 31], this 
expression for tq is twice as large as the corresponding expression from (|4^ , 
although Tq < Hq^ in both cases. For / = 2p/3 however, we find Tq < Hq^ 
from (|46|), while ( pT]) yields Tq > H^^. In the latter case the nonequilibrium 
is 'frozen in', in the former one it is not. This demonstrates explicitly that 
under certain circumstances particle production may improve the conditions 
for bulk viscous inflation. On the other hand, ( |9l| ) makes sense only for 
f/p > 1/2 since for / = p/2 the relaxation time diverges. Especially the 
case / = p/3, dealt with for F = in section 2, is impossible for F = 3-ffo- 

The case F = 3Hq is singled out in different respects as well. For any 
F < 3Hq the particle number density n is decreasing according to (|6^), while 
according to (^) the temperature T increases correspondingly to guarantee 
po = Hq = 0. In the limit F = SHq, n becomes constant and it follows from 
( p6| ) that the temperature is constant as well. Only in this extreme limiting 
case which corresponds to s = (this is only possible since we restricted 
ourselves to a time interval with F = const), the temperature may remain 
constant in a de Sitter phase driven by an effective bulk pressure. For 
F < 3Hq with s 7^ the temperature neccessarily increases, although for 
F 7^ at a lower rate than for F = 0. 

3.3 Entropy production 



Eliminating vr from (|63| ) with the help of ( |36D and (41) the time dependence 
of the entropy per particle is determined by 

nTs = {p + p)[3H-r]+p, (92) 
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or 



nTs 



N-y 



With (|6|) and 



N-y J R^-y ■ 



from (p6D, where the subscript refers to some initial time, one finds 

7a— 1—6 



1 



Po 



noTo V p 



[r^No 



R^, 



pR^ 



(93) 



(94) 



(95) 



In the inflationary phase with p = po and H = Hq, we have 

■ya—l—b 



Po 



" noTo^\RmoJ 
which for radiation (7 = 4/3, a = 1, b = 1/3) reduces to 



(RiN_y 



[3Ho - r] 



4/00 
noTo 



Ha 



Integration of { p7\ ) yields 
noTo 



Ho {t - to) 



1 /■* 



3 7o 



+ s {to) 



(96) 



(97) 



(98) 



The basic difference compared with the case F = is that the change of the 
entropy in a comoving volume, S = nsR^, is no longer determined by the 
change of s alone. With the equations of state for radiation, S is given by 

S = |i^ Ho{t-to)-\f^Tdt +sito)^Noexpj\dt (99) 

in the inflationary phase. If T is again assumed to be approximately given 
by r = SXHq, the expression (99) reduces to 

4po 



E(t) 



noTo 



[1 - A) Ho {t -to) + s (to) 



N (to) exp [3XHo (t - to)] . (100) 



It follows that for F 7^ we have an exponential increase of the comoving 
entropy S during the de Sitter phase. Starting from the latter expression for 
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T, it is possible to apply Maartens' numerical estimations [29] concerning the 
entropy production during the inflationary phase. All his considerations of 
this point are valid in the present case for A 7^ 0, provided his Hq is replaced 
by XHq. His conclusion that it is possible 'to generate the right amount of 
entropy without re-heating' may hold in a universe with particle production. 
There is complete equivalence to Maartens' result for the entropy production 
for A = 1. Consequently, a substantial production of entropy during a (bulk 
viscosity induced) dissipational inflationary phase is possible if at least a 
part of the bulk viscous pressure is related to an increase in the particle 
number. In a sense we have reestablished Maartens' result, although within 
a different setting: An exponential growth of the comoving entropy is only 
possible if the viscous pressure vr in (||) is, at least partially, a 'creation' 
pressure connected with an increase in the number of fluid particles rather 
than an increase in the entropy per particle. 

It remains to consider the stability of the inflationary solutions for F 7^ 0. 
Proceeding as in section 2 we have 

h + (SHoK - TM) h = Q (101) 

instead of (|59[), where 

M = 2—-^ + vl-h. (102) 

IPP + P 

Consequently, the inflationary solutions are stable for 3HqK — TM > 0. 
For n = and the equations of state for radiation one has M = and the 
stability condition reduces to K > again. The stability properties of the 
corresponding inflationary solutions are not affected at all by a nonvanishing 
r. With equations of state for matter the inflationary solutions have the 
same instability as for T = 0. 

4 Summary 

Using particle number density and temperature as basic thermodynamical 
variables of the cosmic fluid we have studied the full Miiller-Israel-Stewart 
theory for a spatially flat, homogeneous and isotropic universe with bulk 
viscous pressure. We found general criteria for the applicability of the so 
called truncated versions. It was shown that in exceptional cases there exist 
common solutions of the full and the truncated theories far from equilibrium. 
The possibility of exponential inflationary solutions driven by bulk viscosity 
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was investigated. Almost all corresponding solutions imply an exponential 
growth of the temperature during the de Sitter stage. Since the number 
density of the fluid particles decreases exponentially, only a corresponding 
increase in the temperature guarantees a constant energy density as long 
as the specific heat is positive and finite and the energy density increases 
with the particle number density. In the second part of the paper the bulk 
viscous pressure was allowed to account partially or fully for particle pro- 
duction processes. A nonvanishing particle production rate may enlarge the 
relaxation time for a viscous pressure that is supposed to drive inflation. It 
may help to 'freeze in' the corresponding nonequilibrium, i.e., to improve the 
conditions for inflation. There exist stable, inflationary solutions for equa- 
tions of state close to that for radiation. Due to the increase in the particle 
number the comoving entropy increases exponentially in this period. Only 
in the limiting case that the entire inflation driving viscous pressure is due 
to particle production both the number density and the temperature remain 
constant during the de Sitter phase. 
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